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Abstract 

We compute the ratio of the inclusive jet cross-sections obtained with the same jet algorithm 
at two different values of the jet radius. We perform a computation of that observable at 
NLO (O {ctg)) in perturbative QCD and compute non-perturbative corrections from soft-gluon 
emission. We discuss predictions for RHIC and the LHC. 

Introduction In this letter, we are interested in computing the ratio of the inclusive jet cross- 
section computed with the same jet algorithm at two different values of R: 

n{pt;Ri,R2) = ^ -. (1) 

Our main objective is to show that the minimal effort required to get a reliable prediction for TZ is 
to include O (a^) perturbative correction as well as (universal) non-perturbative power corrections. 

We shall first discuss the perturbative computation of TZ. This is interesting per se since, as we 
will see below, instead of computing naively the ratio of the cross-sections computed at NLO, which 
would formally correspond to a computation of 7^ up to O (a^), it is actually possible to obtain the 
O (q^) corrections 

Since jets are basic observables at the LHC, and both ATLAS and CMS plan to use the anti-Zct 
algorithm [3] with two different radii (R = 0.4 and 0.6 for ATLAS, R = 0.5 and 0.7 for CMS), they 
could in principle measure the ratio. Compared to the measurement of the inclusive jet cross-section, 
the ratio would not have the uncertainty on the luminosity measurement and would probably be 
less sensitive to the jet energy scale. Below, we shall compare our NLO QCD predictions with and 
without hadronisation corrections to the ATLAS recent measurements |4J and make predictions for 
the cross-section ratio. 

Another situation in which the cross-section ratio is an interesting observable is at RHIC where 
it can be measured both in proton-proton and heavy-ion collisionqj. Due to the interaction with 
the hot medium produced in heavy-ion collisions, one expects the jet to loose energy and to be 
broadened. That would directly translate into a decrease of the cross-section ratio (see e.g. for 
a computation with and without medium effects at O (as) in QCD). Here we shall show that, for 
the pp reference measurement, the next order and hadronisation bring large corrections. 



^See also [2] for an experimental measurement and QCD computation of jet cross-section ratios with different jet 
algorithms and a fixed R. 

^Similar considerations would hold for PbPb collisions at the LHC, with the extra complication that the energy 
of the collision differs from the pp one. 
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Perturbative expansion Let us start witli a perturbative QCD computation of the cross-section 
ratio. Naively, since inclusive jet cross-sections are known up to NLO accuracy (O (af)), one would 
us^ 

n{pt,R^,R2)- ^^^0(^p^.R^y (2) 

which is formally an O (os) computation of TZ. 

The interesting point is that, by making the perturbative expansion explicit, TZ can actually be 
computed up to O (a^). To see this, consider the perturbative expansion of the jet cross-section: 

a{pt; R) = ay^\pt) + a^^^a^pt; R) + afa^^) {pf, R) + (a^) , 

where we have taken into account the fact that the leading-order contribution does not depend 
on R. The contribution a'^'^\pt;R) at a given order a" receives contributions from tree diagrams 
with n particles in the final-state, up to (n — 2)-loop diagrams with 2 particles in the final state. 
Denoting by a^'^'P\pt] R) the p-loop contribution to a^"'\pt; R), we have 

a{pt;R) = ay^'^Hpt) + al[a^'''Hpt;R)+a^'''Hpt)] 

+ at[a(^'^\pt; R) + a^^'^Hpf, R) + a^^'^Hpt)] + O (a^) , (3) 

where we have again used the fact that the contributions with only 2 particles in the final state do 
not depend on the jet radius R. 

If we use to expand the ratio 7^(i?i, i?2) in series of ag, we obtain, up to corrections of order 

-vt T,^ 1^ Aa(^'^)ipt;Ri,R2) 

TZ{pt;Ri,R2) = l + os 7^7—- (4) 

2 Aa(4-Q)(pt; Ri,R2) + A(7(3.i) (p,- R^,R^) y^)ipt;R2)Aa(-''^\pt;Ri,R2) 
a^ipt) [ct(2)(p,)]2 

with Aa("'P)(pt;i?i,i?2) = -a("'P)(pt;i?2). 

The remarkable fact, that allows for the computation to be performed at O (a^) , is that the 
two-loop contribution to the NNLO jet-cross-section, that prevents one from obtaining a NNLO 
computation of the inclusive jet cross-section {e.g. using NLOJet [6J), does not appear in the 
computation of the cross-section raticlf|. 

In what follows, the LO ratio will refer to (jH with the two first terms kept — the O {as) 
expansion, i.e. the first non-trivial order — , while the NLO ratio will also incorporate the O (a^) 
corrections in dH). 

Before proceeding with the discussion about non-perturbative effects, it is interesting to com- 
ment a bit on eq. In the collinear limit, the NLO (resp. NNLO) correction to the cross-section 
will be proportional to Qslog(l/i?) (resp. log^(l/iZ)), which would be the dominant correction 
at small jet radius. In the computation of the ratio, the O (a^) term only involves the cross- 
section difference and will thus be proportional to log(i?i/i22) while the next order will involve 
log(i?i/i?2) log(l/i?2)- This means that for i?i ~ i22 ^ 1, the collinear contribution will mostly 
appear from NLO onwards and we may thus expect large NLO corrections. 



•^For readability, we use a{pt; R) as a shorthand notation for the difTerential cross-section ^{R). 
*Note that it would contribute at the next order. 
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Non-perturbative corrections As we shall see later when making explicit computations of 
the ratio TZ, for small values of R, hadronisation corrections may have a significant impact on 
the jet cross-section and thus on TZ. One could in principle rely on Pythia [7] or Herwig [8] (or, 
better, a combination of both) in order to estimate the correction factor one has to apply to go 
from a parton- level cross-section to a hadron-level cross-section, i.e. to estimate hadronisation 
corrections. Keeping in mind that we want to provide as simple a description of the cross-section 
ratio as we can, we shall instead give an analytic estimate of the hadronisation corrections. In [1], 
hadronisation correction are computed from soft-gluon emission and the authors obtain that the 
effect of hadronisation is to shift the pt of the jet by an average amount 

('^Pt/hadr = — ^ —Aifll). (5) 

In that expression, if the Casimir factor which should be for quark jets and Ca for gluon 
jets, M is the Milan factor that depends on the jet algorithm — it is universal [9], M w 1.49, for 
the anti-kt algorithm while, for the kt algorithm, one finds [lUI M « 1.01 — , and A{fii) carries all 
the non-perturbative dependence. The latter can be rewritterp as 



IT 



ao(/i/) - asipt) - — l^log l^—j + — + l]a^{pt 



(6) 

where the average coupling in the infrared region ao(/i/) = {l/^i) Jq' as{kt)dkt is frequently en- 



countered in event-shape studies (see e.g. fH]), /3o = (HCa — 2n/)/3 and K = Ca y% — "f'^/- 
Including the hadronisation corrections to the perturbative cross-section can then be done usin^ 

(^.r.^ (y{Pt;R) _ ^pQCD(Pt - ('^Pf)hadr;-^) _ ^Lo(pt - {5pt)^^A,]R) 
J^hadr[Pt,R) = 7 57 ~ 7 ~ 7 ^ • [1) 

o"pQCD {pt;R) 0-pQCD {pt;R) (^lo [Pt ; R) 

For the first equality, we have neglected the dispersion in 5pt {i.e. assumed that the shift was always 
the average one) which would correspond to higher power corrections that are not as well controlled 
from LEP data. Approximating the full perturbative cross-section by the leading-order expression 
in the second equality is motivated by the fact that the computation of hadronisation corrections 
from soft-gluon emission is done for the underlying 2—7-2 scattering i.e. from the leading-order 
process. 

Finally, the cross-section ratio after taking into account the hadronisation corrections is 

n{pt;RuR2) = ^^^^f^,n^QCB{Pt;RuR2), (8) 

-f^^hadr (PU R2) 

with TlpQCB {pt'i Rii R2) computed from eq. 

Because of the 1 /R behaviour of ([5]) , we may also expect sizeable effects from the non-perturbative 
corrections at small R. Note however that the factor of the 1/R term is rather small (2Cir^(/x/) w 
0.5 GeV), compared to the corresponding QCD corrections that would typically scale like asPt and 
so dominate at moderate R and pt- 



^At the 2-loop accuracy and in the MS scheme. 

^In practice, since quark and gluon-jets have a difTerent pt shift due to hadronisation, one should consider their 
contributions separately. 
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Figure 1: QCD predictions for the ratio 7^(p^; 0.2, 0.4) at RHIC (^s = 200 GeV) for the anti-A;* 
(left) and kt (right) jet algorithms. On the top panel, the solid lines correspond, from top to 
bottom, to the LO QCD computation (green), to the NLO QCD ratio (red) and to the NLO QCD 
computation including hadronisation effects (blue). The uncertainties due to the scale choice and, 
when relevant, hadronisation are shown as shaded bands on the top panel and the relative scale 
uncertainty is plotted on the bottom panel. For comparison, we have also plotted in dashed lines 
the parton-level (red) and hadron-level (blue) predictions from Pythia. 



Comparison with experiments In the following lines, we briefly discuss the perturbative com- 
putation of TZ and the hadronisation corrections at two different energies: RHIC {^/s = 200 GeV) 
and the LHC (^i = 7TeV). 

As far as the perturbative part of the computation is concerned, we have used NLOJet (v4.1.2) 
[6] for the computation of the different pieces in We have considered the CTEQ6.6 NLO PDF 
set [H] as well as the MSTW08 NLO and NNLO sets pLSj though, for brevity, we shall only show 
the CTEQ6.6 results in what follows. The scale uncertainties have been obtainec|l| by varying 
independently the renormalisation and factorisation scales from pt,jet to ptjot/2 and ^Pt, jet- 

To compute the hadronisation corrections, the only parameter we neeqfl is ao- As already 
mentioned, this can be extracted [llj from event-shape distributions at LEP and we shall use the 
value ao(;U/) = 0.503 with fij = 2 GeV, obtained from JADE data [14] . The uncertainty on the 
hadronisation corrections will be estimated by varying the Milan factor (M = 1.49 for the anti-Zct 
algorithm and M = 1.01 for kt) by the standard 20%. 

Let us start by discussing the case of RHIC, where STAR is planning to measure [15j the ratio 
7^(pt; 0.2, 0.4) for both the kt [16J and anti-fet [3] algorithms. Though the ratio will be measured in 
proton-proton and gold-gold collisions with the ultimate goal to see jet-broadening effects due to 

'^For both scales we compute a negative and a positive uncertainty. The renormalisation and factorisation scale 
uncertainties are then added in quadrature to obtain the total uncertainty. 

*We will always consider large-enough pt so we can safely use n/ = 5 in © and, for consistency, we have used the 
running coupling provided together with the PDF set. 
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Figure 2: Left: comparison of our inclusive jet cross-section computations with the ATLAS mea- 
surements. The ratio between the experimental values and the theory predictions (NLO QCD 
including non-perturbative effects) is plotted; the yellow band represents the uncertainty on the 
theoretical prediction and the cyan band corresponds to the NLO QCD prediction without hadro- 
nisation corrections. Right: predictions for the ratio 7^(pt; 0.4, 0.6) for the anti-fct algorithm at 
,from top to bottom, LO, NLO, NLO with hadronisation corrections. See Fig. [T]for conventions. 



interaction with the hot medium produced in heavy-ion collisions, we just focus on the pp case herll. 
The result is presented in Fig. [1] for both algorithms. The first message is that NLO corrections 
to TZ are substantial (~ 0.1 — 0.15) and, probably as a consequence, the scale uncertainty does not 
decreases when going from LO to NLO. Though they are strictly the same only at LO, the kt and 
anti-Zct algorithms show a very similar cross-section ratio also at NLO. Then, as a consequence of the 
choice of rather small values of R, hadronisation effects are also sizeable (~ 0.15 — 0.2). In this case, 
since the Milan factor is a bit larger for the anti-fct algorithm than for kt, the final ratio tends to be 
a bit larger for the kt algorithm. Finally, Fig. [1] shows that the NLO pQCD computation of 7^ is in 
good agreement with what is obtained from Pythig0 (v6.4) at parton- level (i.e. including parton 
shower from initial and final-state radiation), and our final prediction, including non-perturbative 
corrections is also in good agreement with what Pythia predicts when hadronisation is includec0. 

We now turn to the case of measurements at the LHC and, more precisely, to the jet cross-section 
measured very recently at ^/s = 7 TeV by the ATLAS collaboration [1] . On Fig. [2] we have plotted 
our predictions both for the jet cross-section (anti-kt algorithm with R = 0.4) and the cross-section 
ratio TZ{pt] Ri = 0.4, = 0.6). The jet cross-section is compared to the ATLAS measurements 
and we see that, though the pure NLO QCD prediction (cyan band) describes the data nicely, the 
inclusion of the non-perturbative power corrections improves the description. Note also that the 



^See [5] for a LO description of TZ for pp and gold-gold collisions, incorporating medium effects for the latter. 
^°For Pythia simulations, the ratio is obtained by explicitly dividing the jet cross-section computed with the two 
radii. 

^^The underlying-event corrections could also be taken into account both in our computation and in the Pythia 
simulation but they have a very small impact on TZ. 
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non-perturbative corrections obtained in our approach are compatible with the numbers obtained 
from Pythia and Herwig and quoted by ATLAS. If we now consider the cross-section ratio, see the 
right plot on Fig. [21 we basically recover the main features already discussed in the case of RHIC. 
However, both the NLO QCD corrections and the hadronisation corrections are reduced compared 
to what we observed at RHIC. This is even more true for the non-perturbative corrections at 
large pt which become very small. This is likely due to two effects: first, the considered radii are 
larger, reducing the effect of the collinear divergence in the NLO QCD computation as well as the 
hadronisation corrections that behave like 1/R. Then, the inclusive jet cross-section is much less 
steep at the LHC than at RHIC and thus a common pt shift would have a larger impact at RHIC. 

Conclusions To summarise, we have discussed in this letter the minimal ingredients needed to 
get a reliable calculation of the ratio Tl{pt;Ri,R2) of the pt-dependent inclusive jet cross-section 
computed with the same jet algorithm at two different values, Ri and R2, of the jet radius. 

We have seen that by making an explicitly expansion in powers of a^, we can compute TZ 
perturbatively at O (a^), the NLO accuracy for that observable, that is one order higher than what 
we would naively expect from the direct ratio of the cross-sections. The explicit computation of TZ 
at NLO can be done e.g. using the NLOJet-|-+ event generator. Note that using techniques of |17] 
would allow us to obtain an approximate NNLO calculation and further test the convergence of the 
perturbative series. 

Then, we have estimated the non-perturbative corrections to the ratio. They are based on 
universal power corrections and the only free parameter, qq, can be estimated from fits to event- 
shape measurements at LEP. 

Finally, we have seen that, in practice, both the O (a^) terms and the non-perturbative effects 
are numerically sizeable, except for the hadronisation correction at large pt. In the case of the 
recent jet measurements done by ATLAS, it would be interesting to see if the computation of the 
ratio could benefit from reduced uncertainties compared to the jet cross-section itself. 

In the case of RHIC, the NLO and hadronisation corrections are even larger. It is important to 
keep that in mind when performing the same computation for heavy-ion collisions, in the presence 
of the medium: the one-gluon-emission approximation is likely to be insufficient. One has to include 
the next-order corrections as well as non-perturbative effects. 
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